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PREFACE 


The  attitude  of  a  satellite  in  orbit  is  materially  affected  by 
the  gravitational  torques  which  are  exerted  upon  it.  The  response  be¬ 
havior  of  this  vehicle  to  these  torques  will  depend  on  both  the  initial 
orientation  and  angular  velocity,  and  the  orbital  path  followed. 

Past  studies,  in  which  the  distributed  mass  of  the  satellite  has 
been  approximated  by  a  dumbbell-shaped  body  contained  in  the  orbital 
plane,  and  moving  along  a  circular  orbit,  have  revealed  the  essential 
features  of  the  attitude  response  behavior. 

This  Memorandum  extends  the  analysis  of  the  dumbbell's  motion  to 
elliptic  orbits  and  shows  to  what  extent  the  conclusions  reached  in 
the  circular  case  are  affected  when  orbital  eccentricity  is  taken  into 
account . 
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SUMMARY 

This  Memorandum  discusses  and  analyzes  the  first-order  effect  of 

orbital  eccentricity  on  the  planar  tumbling  or  oscillatory  motion  of 

a  dumbbell -shaped  satellite.  This  has  been  done  by  assuming  that  the 

angular  orientation  angle  Y  can  be  represented  by  a  power  series  in 

eccentricity  e,  in  which  the  coefficient  of  the  e°  term  was  set  equal 

to  the  circular  solution  ¥  available  from  earlier  investigations.  The 

c  1 

differential  equation  for  the  coefficient  of  the  e  term  is  shown  to 
be  of  an  inhomogeneous  Mathieu  type,  the  particular  solutions  of  which 
can  be  readily  obtained  if  certain  weak  restrictions  are  placed  on  the 
initial  magnitude  of  the  dumbbell's  angular  velocity. 

The  analysis  indicated  that  the  orientation  of  the  satellite  in 
the  elliptic  orbit  can  differ  substantially  from  the  one  determined 
for  the  circular  orbit. 
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SYMBOLS 

coefficients  of  the  series  expansion  of  the  Mathieu 
functions 

semimajor  axis  of  the  orbit 
center  of  gravity 
center  of  mass 

Mathieu  functions  of  fractional  order 

eccentricity 
force 

net  force  exerted  . . 

quarter  period  of  the  elliptic  function 

modulus  of  the  elliptic  function 

radial  distance 

half  length  of  dumbbell 

distributed  mass 
total  mass 

mean  anomaly,  n(t  -  tQ) 
mass 

mean  angular  velocity  in  the  elliptic  orbit, 

coefficient  of  trigonometric  term  in  Mathieu's  equation 
parameter  defined  below  Eq.  (43b) 
orbital  radius 

independent  variable,  s  =  e»M 

kinetic  energy 

time 

time  of  perigee  passage 
potential  energy 


X 


V  ■  velocity. 

W  ■  Wronskian 

z  =  independent  variable  defined  by  Eq.  (44) 

a  =  average  value  of  angular  velocity  dY  /dM(>3) 

c 

0  =  order  of  the  Mathieu  function 

^  =  earth's  gravitational  constant,  GMe 

T  -  orbital  period 

9  =  polar  position  angle  of  dumbbell's  C.M. 

Y  -  amplitude  of  oscillation 

Y£  *  amplitude  of  oscillation  in  circular  orbit 

Y  =  initial  angular  velocity  of  dumbbell 
Ci 

U)  *  absolute  angular  velocity  of  dumbbell 
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I.  INTRODUCTION 


An  object  of  arbitrary  mass  distribution  whose  center  of  mass 

2 

moves  in  a  given  orbit  through  a  central  inverse  r  force  field  will 
experience  a  disturbing  torque  acting  around  its  center  of  mass,  which 
would  cause  its  spatial  orientation  to  vary  with  time..  This  torque 
vanishes  when  the  three'  principal  moments  of' inertia  of  the  body  are 
equal.  The  motion  of  a  dumbbell-shaped  object,  being  the  simplest  one  . 
to ' analyze  mathematically,  while  at  the  same  time  exhibiting  many  of 
the  characteristics  which  result  from  the  action  of  this  imbalancing 
gravitational  torque,  has  received  the  most  attention. 

An  understanding  of  the  rotational -motion  characteristics  of  a 
dumbbell -shaped  satellite  could  conceivably  be  of  more  than  purely 
academic  interest,  judging  by  the  frequent  mention  of  such  a  concept 
in  connection  with  a  manned  space  station  in  which  gravity  is  artifi- 
cally  simulated.  The  most  recent  of  these  proposals  is  General  Elec¬ 
tric  Company's  Pseudo  Gravity  Space  Station  concept,  which  envisages 
two  inhabited  compartments  hinged  by- means  of  100-ft-long  extension- 
tubes  and  rotating  around  a  centrally  located  nonrotating  hub. 

In  one  of  the  earlier  papers  on  the  subject,  Klemperer  and  Baker, 
drawing. ori  the  long-known  stable  oscillatory  librational  motion  of  the 
moon  around  its  axis  of  rotation  and  being  interested  in  the  possibil¬ 
ity  of.  utilizing  this  effect  for  measurements  of  orientation  on  board 
a  satellite,  analyzed  the  librational  motion  of  a  symmetric  dumbbell 
placed  in  a ' circular  orbit.  Assuming  small  deviations  from  the  equi¬ 
librium  position,  they  show  that  for  a  dumbbell  of  small  linear  dimen¬ 
sions  (compared  with  the  orbital  radius)  the  dumbbell  behaved  essen¬ 
tially  like  ah  undamped  simple  pendulum  with  a  natural  frequency  equal 
to  ,/3  times  the  .frequency  of  the  motion  in  the  circular  orbit.  They 
also  briefly  mention  that  for  an  orbit  which  is  not  precisely  circular, 

the  situation  would  become  more  complex. 

(2) 

Schindler,  confining  his  analysis  also  to  circular  satellite 
orbits  but  lifting  the  restriction  on  the  size  of  the  angular  devia¬ 
tion,  shows  that  the  librational  periods  can  range  in  a  continuous 
manner  anywhere  from  around  1/3  t  (the  orbital  period)  on  up  to  infinite 
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values  of  the  time,  as  the  initial  angle  which  the  dumbbell  makes  with 
the  vertical  ranges  from  near  zero  to  rt/2  and  that,  therefore,  an  ini¬ 
tial  value  of  this  angle  exists  for  which  the  librational  period  just 
equals  the  orbital  period. 

(3) 

Klemperer,  in  an  ARS  Journal  note,  presents  a  first  integral 
of  the  motion  of  a  librating  dumbbell  and  indicates  that  the  maximum 
value  of  the  amplitude  of  oscillation,  ¥,  of  the  dumbbell  is  defined 
by  the  relation  sin  '4'max  =  YQ/^3n,  where  is  the  angular  velocity 
when  ¥  =  0,  and  n  is  the  angular  velocity,  of  the  center  of  mass  in  the 
circular  orbit.  He  points  out,  furthermore,  that  the  motion  remains 
oscillatory  as  long  as  < y/3n  but  becomes  one  of  tumbling  when  . 
exceeds  this  value. 

(4) 

In  the  same  issue  of  the  ARS  Journal .  Baker  considers  an  orbit 
of  small  eccentricity  and  points  out  that  under  the  restrictive  assump¬ 
tion  of  small  librational  amplitude,  the  motion  is  governed  by  a  forced 
Mathieu  equation. 

Schindler^  employs  a  Hamiltonian  formulation  to  study  the  dy¬ 
namic  behavior  of  the  dumbbell  in  the.  vicinity  of  the  two  positions 
of  equilibrium  and  proves  in  a  rigorous  manner  not  only  that  the  ver¬ 
tical  -dumbbe 11  orientation  is  one  of  stable  equilibrium,  while  the 
horizontal  one  is  a  position  of  unstable  equilibrium,  but  also  that 
a  purely  oscillatory  motion  around  the  unstable-equilibrium  position 
is  possible  in  an  elliptic  orbit,  provided  the  initial  conditions  are 
properly  chosen. 

Moran^  employs  a  perturbation  technique  to  study  the  effect 
which  the  librational  mode  of  motion  of  a  dumbbell  moving  initially 
along  an  undisturbed  circular  orbit  has  on  the  motion  of  its  center 
of  mass.  He  shows  that  resonance  conditions  in  the  radial  equation 
of  motion  can  arise  for  certain  values  of  the  initial  librational  fre¬ 
quency  and  that  the  differential  equation  for  the  first  radial  cor¬ 
rection  term  gives  rise  to  a  divergent  oscillatory  solution  as  time 
increases. 

The  present  Memorandum  describes  in  a  qualitative  and  quantita¬ 
tive  manner  the  first-order  effects  of  orbital  eccentricity  on  the 
planar  librational  motion  of  a  dumbbell  by  means  of  a  perturbation  ex¬ 
pansion  in  the  vicinity  of  the  nominal  circular-orbit  solution. 
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II.  THE  GRAVITATIONAL  TORQUE 

Before  setting  out  to  write  down  in  a  mechanical  fashion  Lagrange's 
equations  governing  the  planar  motion  of  the  dumbbell  satellite,  it 
might  be  useful  and  instructive  first  to  consider  briefly  one  aspect 
of  the  physical  origin  and  nature  of  the  gravitational  torque  which'  is 
but  seldom  recognized  and  mentioned.  Reference  to  Fig.  1  makes  it  im¬ 
mediately  obvious  that  as.  mass  m^  of  the  symmetric  dumbbell  shown 
(m^  =  m^  =  m)  is  closer  to  the  center  of  attraction  than  is  mass  m2 , 
it  will  experience  a.  larger  gravitational  radial  force  than  will  mass  1 
nu,. .  The  difference  in  magnitude  and  orientation  of  those  two  forces  .  \ 
gives  rise  .to  the  well-known  expressions  for  the  resultant  force,  which 
acts  through  the  center  of  mass  0  and.  perturbs  the  translational  mo¬ 
tion,  and  to  the  gravitational  torque  around  0,  which  gives  rise  to 
•the  librational  motion. •  .  • 


Fig.  I  —  Geometry  and  notation 


4 


A  more  graphic  way  of  interpreting  the  action  of  the  gravitational 
torque  and  one  which  allows  for  a  more  qualitative  feel  for  its  sense 
and  magnitude  is  by  a  visualization  of  the  instantaneous  relative  po¬ 
sition  of  the  dumbbell's  center  of  mass  (C.M.)  and  center  of  gravity 
'  (C.G.) .. 

For  a  body  of  arbitrary  shape,  these  two  points  coincide  in  a  con¬ 
stant,  unidirectional  force  field,  and  for  spherically  symmetric  bodies, 
2 

.  also  in  an  inverse  r  field.  In  general,  though,  the  two  points  are 
distinct; 

The  center  of  gravity  of  a  distributed  mass  M  is  found  by  an  ap¬ 
plication  of-  Newton's  Third  Law.  It  is  that  point  at  which  a  concen¬ 
trated  mass  equal  to  the  total  mass  (M  =  2m)  would  experience  a  force- 
equal  in  magnitude  and  opposite  in  direction  to  that  experienced  by 
the  earth  at  its  center,  due  to  the  gravitational  attraction  exerted 
upon- it  by  the  mass  M. 

This  definition  is  now  used  to  locate  the  center  of  gravity  of 
the  dumbbell  for  two  extreme  cases  of  its  orientation,  f  =  n/2  and 
Y  =  0,  as  shown  in  Fig.  2. 

4  ®-C.G'.  •  m2 


)C.M. 

C.G. 

m,  r-o 

.'Center  of  earth  • 

(a)  Center  of  earth 

(b) 

Fig,  2 — Relative  location  of  center  of  mass  and  center  of  gravity 
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Letting  F^  and  F£  designate  the  forces  at  and  respectively 
we  have 


m  - 

3  rl 


_  31  “ 

3  r2 
r2 


(1) 


where  =  n^  =  m  and  p,  *•  GMg,  the  earth's  gravitational  constant.  Be 
cause  of  syranetry,  the 
parallel  to  r,  so  that 


cause  of  syranetry,  the  net  force  exerted  on  the  dumbbell,  Fr,  will  be 


F  m  J  •  -7  +F  •  I  ■  -  53i.  C  i.  •  T  ^  •  SUl  T  .  .  T  'N 
r  1  r  2  r  2  V  r,  r  J  2  v  r_  r  J 

r,  1  r.  2  ^ 


2mnr 


(2) 
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This  force  acts  on  the  concentrated  mass  2m  located  at  a  radial 
distance  L.  Under  the  assumption  l/r  «  1  we  get 


from  which  we  find  that 


r  +  Ar 


(4) 


and 


Ar 


2  £ 
4  r 


(5) 
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This  places  the  center  of  gravity  above  the  center  of  mass.  In 
the  case  of  a  dumbbell  rotating  in  a  counterclockwise  direction,  any 
slight  increase  in  the  angle  ¥  beyond  tt/2  would  cause  the  center  of 
gravity  to  shift  to  the  left  of  the  vector  r  and  thereby  give  rise  to 
an  accelerating  moment  around  the  center  of  mass.  The  converse  holds 
true  during  the  period  of  time  prior  to  the  reaching  of  the  horizontal 
prientation,  when  the  torque,  has  ai  retarding  effect. 

A  repetition  of  the  above  analysis  for  the  orientation  (b)  of 

Fig.  2  shows  that  the  center  of  gravity  is  now  located  a  distance 
2 

Ar  ■  3/2  i  /r  below  the  center  of  mass.  It  is  immediately  evident 
that  any  slight  deviation  from  the  vertical  position  gives  rise  to  a 
torque  which  opposes  the  motion,  and  thus  labels  the  vertical  position 
one  of  stable  equilibrium. 

The  path  described  by  the  dumbbell's  center  of  gravity  during  one 
complete  period  of  rotation  is  shown  schematically  in  Fig.  3  and  gives 
an  indication  of  the  oscillatory  nature  of  the  gravitational  torque. 


Fig.  3 — Locus  of  center-of-gravity  travel  (not  drawn  to  scale) 
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III.  EQUATIONS  OF  MOTION 


The  kinetic  energy  is 


■  2»Cvi  +  ,0  (6) 


and  the  potential  energy  is 


U 


^  ^ 


(7) 


where 


r  +  jftui 


V„  -  r  -  iu)i 


U)  m  Q  +  Y 

r  -  ri_  +  r0i 


e 


r,  -  |  r 2  +  X2  -  2rX  cos  *] 


[■ 

£r2  +  A2  +  2rX  cos 


1/2 


1/2 


(8) 


Figure  4  shows  the  various  parameters  and  unit  vectors  used  in 
the  above  relations.  After  Eqs.  (8)  are  substituted  into  Eqs.  (6)  and 
(7)  and  the  expression  for  the  Lagrangian  f-U  is  obtained,  one  can 
derive  the  set  of  differential  equations  shown  below. 


¥  -  r0  - 


.  Ul 

2 


t  X  coa  T 

L  ri 


r  +  l  cos  V 
2 


(9) 


2a4t  [r2®  +  i2<®  +  *>]  "  0 


(10) 
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i  u  urjl  siri  ^ 

i  i 

1  *  (8  +  Y)  1  2 

3  "  3 

Lr2  riJ 

(ID 


Equation  (10)  indicates  that  the  combined  angular  momentum  of  the 
motion  is  conserved. 


Fig.  4  —  Geometry  of  motion 


•  The  bracket  on  the  right-hand  side  of  both  Eq.  (9)  and  Eq.  (11) 
can  be  evaluated  with  the  aid  of  the  relations  for  r^  and  r ^  of  Eq.  (8). 
Realizing  that  t/r  «  1,  we  can  write  the  difference  of  the  inverse 
cubes  of  the  radii  as 
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I _ I_  .  _ 1 _ _ i 


3  3  2  3/2 

2  1  ■*♦(*)  ♦<*)-•* ]  *♦.(*> 

2  3/2 

-  2^)  cos  y] 

a  [l  -  3^)  cos  Y  ”  1  -  3^  cos  f  +  0 

\2  6A 

)  ...  -  cos  Y 

'  T 

After  the  right-hand  side  of  Eq.  (9)  is  also  evaluated,  we  come 
up  with  the  following  set  of  differential  equations: 

r  -  rfi2  «  -  “j[l  +  fCji  (3  cos2'?  -  1) 
r  V  y 

+  ...]  (12) 

dF{t2[4  +  0)!  «♦«]}  -  0 

(13) 

0  +  t  +  f  **3  sin  21f.  -  0  ' 
r  . 

(14) 

It  is  observed  that  the  dumbbell  size  appears  only  through  the 

2 

parameter  (A/r)  .  In  order  to  be  able  to  express  the  orbital  parameters 

2  2 

as  power  series  in  (A/r)  it  is  necessary  to  treat  (A/r)  as  a  con¬ 
stant.  This  can  be  done  for  orbits  of  small  enough  eccentricity,  e, 
when  A/r  a  A/a,  where  a  is  the  semimajor  axis  of  the  orbit.  Follow¬ 
ing  Moran^  by  assuming  expansions  of  the  form 

1  Uif 

nm  0 

r-  1 

(15) 

n  *  0 


r 


10 


there  results  the  zeroth  set  of  equations 


Y  +  f  “3  sin  2*  -  -  (18) 

r 

o 

Equations  (16)  and  (17)  define  of  course  the  undisturbed  orbit 
of  the  center  of  mass  and  have  as  their  solutions  the  series 

2 

rQ  -  a[l  -  e  cos  H  -  |  (cos  2M  -  1)  -  0(e3)...]  (19) 

6o  -  M  +  2e  sin  M  +  |  e2  sin  2M  +  0(e3)...  (20) 

where 

M  ■  n(t  -  tQ)  “  mean  anomaly 
tQ  ■  .  time  of  perigee  passage 
n  “  mean  angular  velocity  In  the  elliptic  orbit 

Now  let 

’  ■  *c+  t  Tn‘°  (2l> 

n  ■  1 

where  Y  denotes  the  orientation  angle  of  the  dumbbell  when  the  center 
c 

of  mass  describes  a  circular  orbit.  Substituting  Eqs.  (19),  (20),  and 
(21)  into  Eq.  (18)  and  collecting  like  terms  in  eccentricity  gives,  to 
linear  terms  in  eccentricity 
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Y  +  \  n2  sin  2Y  -  0  (22) 

c  z  c 

Yc  +  eYx  +  |  n2  [l  +  3e  cos  m]  [sin  2(YC  +  eY^]  -  2en2  sin  H  (23) 

Approximating  sin  2(YC  +  eY^)  by  (sin  2YC  +  2eY^  cos  2Yfi)  and 
using  Eq.  (22)  gives  the  following  differential  equation  for  Y^ 

Y,  +  (3n2  cos  2Y  )  Y,  ■  2n2  sin  M  -  %  n2  cos  M  sin  2Y  (24) 

l  cl  z  c 

We  now  replace  time,  t,  by  the  mean  anomaly,  M,  as  the  independent 

2 

variable,  which  eliminates  the  factor  n  in  Eqs.  (22)  and  (24).  Thus 

Y"  +  i  sin  2Y  -  0  (25) 

c  z  c 

YV  +  (3  cos  2Y  )  Y,  -  2  sin  M  -  §•  cos  M  sin  2Y  (26) 

1  C  1  z  c 

where  (  ) *  denotes  differentiation  with  respect  to  the  mean  anomaly  M. 
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IV.  M3TI0N  IN  THE  ELLIPTIC  ORBIT 

We  shall  now  investigate  the  nature  of  the  solutions  of  Eq.  (26). 

These  depend  on  the  form  of  the  solution  of  Eq#  (25) ,  the  differential  ^ 
equation  of  the  circular-orbit  case.  As  shown  by  Schindler  and  Moran, 
the  solution  to  Eq.  (25  depends  on  the  initial  value  of  Y'c  for  “  0, 
and  is  given  by 


sin 


v  ■  eo  i*cj  <y5  (°8ciiiati°n> 


(27a) 


sin  ¥  =  Sn  f  0 

c  c^  o 


I*'  |  >  (tumbling)  (27b) 


The  modulus  k  of  the  elliptic  function  is  given  by 
k  ■  (oscillation) 


(28a) 


k  »  (tumbling) 


(28b) 


Furthermore,  for  two  values  of  the  initial  angular  velocity,  f'  = 

Ci 

1.66675,  1.78794,  the  periods  of  the  oscillatory  and  tumbling  modes, 
respectively,  will  equal  the  orbital  period  of  motion. 

A  plot  of  Eq.  (27),  taken  from  Ref.  6,  is  reproduced  for  conven¬ 
ience  in  Fig.  5. 


SOLUTION  FOR  LARGE  Y '  (TUMBLING  MOTION) 
ci 

The  solution  of  Eq.  (26)  with  Eq.  (27)  taken  into  account  is  fa¬ 
cilitated  if  the  range  of  possible  values  of  is  broken  up  *nto  two 

separate  regions.  For  large  values  of  ¥ '  ,  say  '  2  3,  Fig.  5  shows 

i  i 

that  the  oscillatory  ripple  around  the  constant -slope  line  of  versus 
0q  has  nearly  disappeared,  so  that  it  is  reasonable  to  approximate  sin  Yc 
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Fig.  5 — Time  history  of  rotational  motion  in  a  circular  orbit  (6) 
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3 

For  a  E  3 ,  q  =  -  X  =  “0.1667. 

la 

The  order  of  the  Mathieu  function,  &,  is  obtained  by  equating  to 
each  other  the  two  functions  vq  and  vq,  which  are  given  by  the  expressions 


-q/<@  +  2 Y 


q 


(e  +  2)2  (p  +  4)2 


-  ...  0(qJ) 


(32) 


q  (P  -  2)2  (p  -  2)2  (P  -  4)2 


+  •  ••  0(q  ) 


For  the  present  case,  a  value  of  |3  5s  0.118  is  found*  The  coef- 
ficients  A^  are  solved  from  the  recurrence  relations 


(2r  +  p)2  -  q(A^2  +  A^)  -  0  (-  <  r  <  -)  (33) 


When  terms  beyond  A^,  and  A?  are  neglected,  the  other  coefficients 


can 


T 

i.  of 

aP. 

0 

AP. 

-4 

AP 

-2 

= 

AP 

s 

z 

AP 

A4 

a 

AP 

o 


-0.03716  A1 


P 


(34) 


0.0003651  A' 


P 


and  the  two  solutions  then  become 


w  „  Ce  .  aP  pO. 0005205  cos  3.882s  +  0.04706  cos. 1.882s 

*11  0.118  o  L 

+  cos  0.118s  +  0.03716  cos  2.118s  +  0.0003651  cos  4.118s  +  .../] 


(35) 
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by  sin  oM,  where  M  is  the  mean  anomaly,  and  a  is  the  average  value  of 
the  angular  velocity  d^/dH  which  is  assumed  to  be  larger  than  or  equal 
to  3.  For  convenience  of  analysis  Of  will  be  considered  an  integer.  With 
the  additional  change  of  independent  variable  s  ■  oM,  Eq.  (26)  becomes 


ds 


( 


cos  2s 


,  s  1 
sin  —  - 

“  2a' 


“  cos  —  sin  2s  (29) 
2  a 


which  is  an  inhomogeneous  Mathieu  equation  of  the  form 


+  Ta  -  2q  cos  2x|z  =■  f(x) 
dx  L  J 


(30) 


with 

a  “  0 


and 

q  "  '  ^“2 

2a 

In  view  of  the  assumption  a  i  3,  all  solutions  to  Eq.  (29)  are 
stable  and  oscillatory  in  nature  and  consist  of  Mathieu  functions  of 
fractional  order.  For  purposes  of  illustration,  the  method  of  obtain¬ 
ing  the  solution  for  the  case  a  *  3  will  be  demonstrated  below. 


NUMERICAL -EXAMPLE  SOLUTION 

Employing  the  notation  of  Ref.  7,  the  two  linearly  independent 
solutions  to  the  homogeneous  part  of  the  equation,  and  Y12>  are 
the  functions  Ce^s,  -q)  and  Se^s,  -q)  which  are  given  by  the  series 

Y11  “  Cep(8>  -<I>  "  ("Dr  A2r  cos  (2r  +  P)s 

r  *»  — co 

(31) 

Y12  -  Seg(s,  -q)  -  J  (-l)r  sin  (2r  +  p)s 

r  -  -• 
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Y12  “  Se0  118  =  Ao  [-0*0005205  sin  3.882s  -  0.04706  sin  1.882s 


+  sin  0.118s  +  0.03716  sin  2.118s  +  0.0003651  sin  4.118s  +  ... 

(36) 

The  complementary  solution  to  Eq.  (29)  is  then 


*1  =  Cl%.118  +  C2Se0.118 


(37) 


The  solutions  are  nonperiodic  but  are  bounded  and  tend  neither 
to  0  or  »  as  s  -•  », 

If  we  now  define  the  set  of  fundamental  solutions  Y^2  suc^ 

that 


?H(°) 


f{x(0) 


?{2(0) 


?12(°) 


(38) 


=  0 


the  complete  solution  is  given  by 

**  C^s)  Yu(s)  +  C2(s)  Y12(s)  +  Yx(0)  ?n(s)  +  Y{(0)  ?12(s) 

(39) 


where 


s  h(?)  *19(?) 

C1(s)  -  -J*  - - ~ - -  d5 


0  W[Y11C5> »  Y12(5)] 


s  h(?)  Y„(5) 

C2(s)  =  f  T - - - n  d? 

.  o  w[yu(?),  y12(S)J. 
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MS)  *  sin  ^  cos  ^  sin  2C 

a2  01  2a2  a 

few  ,  ^12^)^]  “  Wronskian  of  Yjj  and  Yj^ 

Since  all  the  initial  conditions  on  the  motion  are  assumed  to 
have  been  satisfied  by  the  circular- orbit  solution,  Yc>  we  have  Y^(0)  = 
Y{(0)  =  0,  and  the  complementary  portion  of  the  solution  consequently 
vanishes.' 

The  Wronskian  W  will  be  a  constant  because  Eq.  (29)  has  no  first 
derivative  of  Y^,  and  can  therefore  be  pulled  out  from  the  integral 
sign.  Solution  (39)  can  then  be  written 

S’  S  ■ 

W(0)Y1  ■  Se0.118^  I  h^)  Ce0.118^d^  “  Ce0.118^  I  Se0.118^d^ 

o  o 

(40) 

When  all  the  operations  are  carried  out  we  finally  end  up  with 
the  solution 

Yx  -  1.00028  -  1.0565  sin  0.3333s  +  0.10892  sin  1.6667s 

•  -  0.01276  sin  2.3333s  -  0.00124  sin  4.3333s 

-  0.00112  sin  3.6667s  +  _  (41) 

The  result  is  plotted  in  Fig.  6,  which  shows  the  behavior  of  Yc*  Y^, 
and  Y^  as  a  function  of  M  for  one  orbital  period  of  the  dumbbell. 

The  magnitude  of  the  perturbation  caused  by  an  orbital  eccentri¬ 
city  of  e  ■  0.1  is  shown  by  the  dashed  curve.  It  is  apparent  that 
orientation  errors  as  large  as  80  deg  could  occur  during  the  first 
orbital  period. 
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SOLUTION  FOR  SMALL  V  (OSCILLATORY  MOTION) 

As  'i,/  is  reduced  from  a  value  of  about  3  and  approaches  73 » 

the  accuracy  of  the  approximation  ¥  rj?'  .M  becomes  progressively 

c  ci 

worse,  and  the  above  solution  ceases  to  describe  the  motion. 

Neglecting  the  region  around  the  condition  f'  =  ,V3,  which  is 

ci 

only  of  limited  interest,  we  consider  now  the  form  of  the  solution 
of  Eq.  (26)  for  ¥'  <  J\. 

ci 

The  functions  cos  2¥fi  and  sin  2¥c  appearing  in  the  differential 
Eq.  (26)  can  be  expressed  as  functions  of  M  by  means  of  Eq.  (27a). 

cos  2¥  =1-2  sin2¥  = 

c  c 

sin  2¥  =  2  sin  ¥  cos  ¥ 

c  c  c 


1-2  k2sn2y3  M 

■  2k  sn  ,/3  M  dn^3  H 


-2k 


(42) 

d - 73  H 


d(  M) 


In  order  to  bring  Eq.  (26)  into  a  more  tractable  form,  the  func¬ 
tions  sn^  M  and  cn^S  H  are  expanded  in  Fourier  series. 


kK  r..  .  V3  nM  .  ,  3  ^3  ttM  .  ^ _ 

2^  snv3  M  =  — :  sin  - gr  +  •— 3-  sin  2K  +  .  -5  8in  2K 

1-q  1  -  q  1-q 


5-73  nM  . 


l 


-r+1/2 

_3 _ 

r  =  0  1  q 


sin 


(2r  +  1)  73 


nM 


2K 


(43a) 


cn73  M 


q 

K 

K' 

,/2 


I  ~ 

r  =  01  + 
-nK'/K 


_r+l/2 


_2rfl 

q 


COS 


(2r  +  1)  nM 


2K 


K(k)  -  J 

K(k') 

1  -k2 


n/2 


d0 


1  -  k2sin20 


(43b) 
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For  sufficiently  small  values  of  q,  the  retention  of  the  leading  tern 
in  each  of  Eqs.  (43a)  and  (43b)  might  give  sufficient  accuracy.  A 
rough  idea  of  the  size  of  q  is  obtained  by  taking  for  instance  the 
case  Y7  “  1.66  as  a  rough  indication  of  an  upper  bound  on  the  angular 

C4 

velocity.  This  corresponds  to  equality  in  orbital  and  librational 
periods.  For  this  case 


-  0.958 

s? 

V-/ 

1 

fo 

• 

00 

0.2868 

K(k')  a:  1.6 

_  1.6 

e  2,8  =  0.165 

As  Y 7  decreases,  k  will  decrease  and  k'  will  increase,  so  that 

cj 

K  will  be  reduced  and  K'  increased;  this  will  lead  to  a  decrease  in 
the  size  of  q.  Thus,  for 


Y'  <  1.66 
Ci 

q  <  0.165 


Placing  an  upper  limit  of  1.66  on  Y  ,  it  appears  that  the  reten- 

Ci 

tion  of  only  one  term  in  the  Fourier  expansions  would  b"  satisfactory 
in  the  present  case.  Under  this  approximation  terms  which  are  of  mag¬ 
nitude  q  have  been  neglected  in  comparison  with  unity.  Thus 


cos  2YC  r-s  1  -  P2  sin2 


.  P2  .  P2  V3  ttM 
1  -  r  +  cos  ~ 


2tt 


2  Vo 


'C  kK2  (1  +  q) 


^3  nH 
2K 


»IT2qY^2 

P2-  - i 


8n2q 


3k2K2(l  -  q)2  K2(l  -  q)2 
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(Note  that  k 


by  definition.) 


Making  the  additional  transformation  of  independent  variables 


z  = 


V3  ttM 


2K 


(44) 


we  obtain  Eq.  (26)  in  its  final  form 


2 

d  Y 


dz 


—  +  j^a  -  2q  cos  2zJ  Y^  = 


81 C  2K  _  6  ^ 


sin 


3tt 


z  - 

•v/3  tt  k(l  +  q) 


^-rsinr-^-+i 

+  a)  L 


+  sin 


( 


2K 


V3  TT 


>] 


\/T" 

(45) 


where 


4Kl  _  q 

2  16  ,  -.2 

tt  (1  -  q) 


(1  -  q)‘ 


As  in  the  last  case,  here  too  the  motion  of  libration  is  governed 
by  an  inhomogeneous  Mathieu  equation.  The  bounded  oscillatory  form  of 
the  forcing  function  suggests  immediately  that  for  finite  values  of 
Y^  the  nature  of  the  forced  solution  of  Eq.  (45)  depends  on  the  sta¬ 
bility  or  instability  exhibited  by  the  solution  to  the  homogeneous 
Mathieu  equation.  This  question  is  easily  resolved  with  the  aid  of 
Fig.  7,  which  presents  in  the  plane  of  (a)  versus  (q)  the  stable  and 
unstable  regions  of  the  solution  to  Mathieu1 s  equation.  Due  to  sym¬ 
metry  with  respect  to  the  a  axis,  the  right  side  of  the  Mathieu  plane 
has  been  used. 

The  curves  CeQ,  Se^ ,  and  Ce^  denote  respectively  the  loci  of  zeroth 
and  first-order  periodic  solutions  to  Mathieu' s  equations  and  represent 
the  transition  curves  between  the  stable  and  unstable  regions  shown. 
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The  directed  dashed  curve  indicates  the  trace  of  the  (a,q)  point  of 
Eq.  (45)  as  the  initial  value  of  angular  velocity  Y^  tends  toward 
zero.  Except  for  a  narrow  region  around  a  =  1,  q  =  & ,  the  homogeneous 
solution  is  at  all  times  oscillatory  and  divergent  as  z  increases, 
making  the  complete  solution  unstable.  That  segment  of  the  dashed 
curve  which  appears  to  follow  closely  along  the  Se^  curve  has  associ¬ 
ated  with  it  divergent,  though  periodic,  homogeneous  solutions. 

Although  the  scale  of  Fig.  7  does  not  make  this  apparent,  a  simpli 
fied  analysis  shows  that  the  dashed  curve  approaches  the  point  a  =  1, 
q  **  0  with  a  slope  da/d|q|  =  -2.  Since  the  Se^  curve  exhibits  at  this 
point  a  slope  of  -1,  it  appears  that  the  dashed  curve  penetrates  into 

the  stable  region  of  the  plane  as  Y'  assumes  vanishingly  small  values. 

Ci 

In  the  limit,  as 

Y '  -  0 

ci 

K  -»  tt/2 
a  -»  1 

—  -  4 

■A 

Equation  (45)  is  reduced  to  the  form 


(46) 


The  complementary  solution  of  Eq.  (46)  will  be  bounded  and  perio¬ 
dic,  but  the  forced  solution,  though  bounded  too,  will  not  exhibit  any 
periodicity.  The  form  of  the  arguments  of  the  trigonometric  terms  in 
the  forcing  function  precludes  the  occurrence  of  resonance  conditions. 
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Fig.7 — Effect  of  initial  angular  velocity 
^  on  the  location  of  the  solution  point  in 
1  the  a,  q  plane 
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VI.  CONCLUSION 


The  planar  rotational  motion  of  a  dumbbell-shaped  satellite  whose 
center  of  mass  moves  along  an  undisturbed  elliptic  orbit  of  slight  ec¬ 
centricity  has  been  analyzed.  It  has  been  shown  that  the  correction 
terms  due  to  eccentricity,  which  must  be  superimposed  on  the  angular- 
position  time  behavior  of  the  same  satellite  when  moving  along  a  cir¬ 
cular  orbit,  are  always  bounded  and  oscillatory  for  all  values  of  ini¬ 
tial  angular  velocity  in  excess  of  roughly  3  times  the  orbital  angular 
velocity.  Depending  on  the  magnitude  of  the  eccentricity,  these  per¬ 
turbations  could  be  quite  appreciable. 

When  the  initial  motion  is  one  of  pendulous  oscillations,  the 
first-order  effect  of  eccentricity  is  to  introduce  a  divergent  oscil¬ 
latory  term  into  the  time  behavior  of  the  first  perturbation  term  and 
thus  to  cause  the  orientation  angle  in  the  elliptic  orbit  to  differ 
significantly  from  that  assumed  in  the  circular  orbit.  Because  of 
possible  phase  differences  in  the  time  behavior  of  the  circular  orien¬ 
tation  angle  Y and  the  perturbed  correction  term  Y^,  it  is  still  possi¬ 
ble  for  the  complete  solution  not  to  exhibit  an  actual  boundless  in¬ 
crease  in  amplitude  as  long  as  the  restrictions  imposed  by  linearity 
are  not  violated. 


